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equilibrium many-boson systems embedded in a thermal bath, better referred-to 
as Frohlich-Bose-Einstein condensation. The complex behavior emerges after a 
threshold of the exciting intensity is attained. It is inhibited at higher intensities 
when the magnon-magnon interaction drives the magnons to internal thermal- 
ization. The observed behavior of the relaxation to equilibrium after the end 
of the pumping pulse is also accounted for and the different processes fully de- 
scribed. 
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1. Introduction 

The kinetic of evolution of the system of spins in thin films of yttrium-iron- 
garnets in the presence of a constant magnetic field, and being excited by a source 
of rf-radiation which drives the system towards far-removed from equilibrium 
conditions, has been reported in detailed experiments performed by Demokritov 
et al. [BE]- These experimental results have evidenced the occurrence of an 
unexpected large enhancement of the population of the magnons in the state lowest 
in energy in their energy dispersion relation. That is, the energy pumped on the 
system instead of being redistributed among the magnons in such non-thermal 
conditions is transferred to the mode lowest in frequency (with a fraction of course 
being dissipated to the surrounding media). Some theoretical studies along certain 
approaches has been presented by several authors (see for example Refs. [31 HI |5]); 
we proceed here to describe the phenomenon within a complete thermo-statistical 
description ensemble formalism. 

Such phenomenon has been referred-to as a non-equilibrium Bose-Einstein 
condensation, which then would belong to a family of three types of BEC: 

The original one is the BEC in many-boson particle systems in equilibrium at 
very low temperatures, which follows when their de Broglie thermal wave length 
becomes larger than their mean separation distance, and presenting some typical 
hallmarks (spontaneous symmetry breaking, long-range coherence, etc.). Aside 
from the case of superfluidity, BEC was realized in systems consisting of atomic 
alkali gases contained in traps. A nice tutorial review is due to A. J. Leggett [6] 
(see also [7]). 

A second type of BEC is the one of boson-like quasi-particles, that is, those 
associated to elementary excitations in solids (e.g. phonons, excitons, hybrid 
excitations, etc.), when in equilibrium at extremely low temperatures. A well 
studied case is the one of an exciton-polariton system confined in microcavities (a 
near two-dimensional sheet), exhibiting the classic hallmarks of a BEC [8]. 

The third type, the one we are considering here, is the case of boson-like 
quasi-particles (associated to elementary excitations in solids) which are driven 
out of equilibrium by external perturbative sources. D. Snoke [H] has properly 
noticed that the name BEC can be misleading (some authors call it "resonance" , 
e.g. in the case of phonons [10]). and following this author it is better not to 
be haggling about names, and we introduce the nomenclature NEFBEC (short 
for Non-Equilibrium Frohlich-Bose-Einstein Condensation for the reasons stated 
below). As noticed, here we consider the case of magnons (boson-like quasi- 
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particles), demonstrating that NEFBEC of magnons is another example of a 
phenomenon common to many-boson systems embedded in a thermal bath (in 
the conditions that the interaction of both generates non-linear processes) when 
driven sufficiently away from equilibrium by the action of an external pumping 
source and which display possible applications in the technologies of devices and 
medicine. 

1. A first case was evidenced by Herbert Frohlich who considered the 
many boson system consisting of polar vibration (lo phonons) in 
biopolymers under dark excitation (metabolic energy pumping) and 
embedded in a surrounding fluid[HJ [T21 H21 EH- From a Science, 
Technology and Innovation (STI) point of view it was considered 
to have implications in medical diagnosis [15]. More recently has 
been considered to be related to brain functioning and artificial 
intelligence [H]. 

2. A second case is the one of acoustic vibration (ac phonons) in biological 
fluids, involving nonlinear anharmonic interactions and in the presence 
of pumping sonic waves, with eventual STI relevance in supersonic 
treatments and imaging in medicine|17[ ITS] . 

3. A third one is that of excitons (electron- hole pairs in semiconductors) 
interacting with the lattice vibrations and under the action of rf- 
electromagnetic fields; on a STI aspect, the phenomenon has been 
considered for allowing a possible exciton-laser in the THz frequency 
range called "Excitoner" [131 120] . 

4. A fourth one is the case of magnons already referred to [TJ |2J, which we 
here analyze in depth. The thermal bath is constituted by the phonon 
system, with which a nonlinear interaction exists, and the magnons 
are driven arbitrarily out of equilibrium by a source of electromagnetic 
radio frequency [21]. Technological applications are related to the 
construction of sources of coherent microwave radiation 



There exist two other cases of NEBEC (differing from NEFBEC) where the 
phenomenon is associated to the action of the pumping procedure of drifting 
electron excitation, namely, 

5. A fifth one consists in a system of longitudinal acoustic phonons 

driven away from equilibrium by means of drifting electron excitation 
(presence of an electric field producing an electron current), which has 
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been related to the creation of the so-called Saser, an acoustic laser 
device, with applications in computing and imaging [24"] [25] . 

6. A sixth one involving a system of LO-phonons driven away from 

equilibrium by means of drifting electron excitation, which displays a 
condensation in an off-center small region of the Brillouin Zone [26|l27j. 

We describe here item number 4, namely, a system of magnons excited by an 
external pumping source. For that purpose, we consider a system of N localized 
spins in the presence of a constant magnetic field, being pumped by a rf-source 
of radiation driving them out of equilibrium while embedded in a thermal bath 
consisting of the phonon system (the lattice vibrations) in equilibrium with 
an external reservoir at temperature T . The microscopic state of the system 
is characterized by the full Hamiltonian of spins and lattice vibrations after 
going through Holstein-Primakov and Bogoliubov transformations [281 l29| 130] . On 
the other hand, the characterization of the macroscopic state of the magnon 
system is done in terms of the Thermo-Mechanical Statistics based on the 
framework of a Non-Equilibrium Statistical Ensemble Formalism (NESEF for 
short) [3U |32j [33J EU ESI [36] . Other modern approach consists in the use 
of Computational Modeling [37] 138] (developed after Non-equilibrium Molecular 
Dynamics [39J). It may be noticed that NESEF is a systematization and an 
extension of the essential contributions of several renowned authors following the 
brilliant pioneering work of Ludwig Boltzmann. The formalism introduces the 
fundamental properties of historicity and irreversibility in the evolution of the 
nonequilibrium system where dissipative and pumping processes are under way. 

In terms of the dynamics generated by this full Hamiltonian the equations of 
evolution of the macroscopic state of the system are obtained in the framework 
of the NESEF-based nonlinear quantum kinetic theory [SJ E21 E31 E2 133 ESI EH 
HU 22J H3]. We call the attention to the fact that the evolution equations are 
the quantum mechanical equations of motion averaged over the nonequilibrium 
ensemble, with the NESEF-kinetic theory providing a practical way of calculation. 

The paper is organized as follows: 

In section 2 the Theoretical Background is described; In section 3 is presented 
the Evolution of the Nonequilibrium Macrostate of the Magnon System; In 
section 4 NEFBEC in YIG is studied; In section 5 the Decay of the Condensate 
is analyzed; Finally, in section 6 we present Additional Considerations and 
Concluding Remarks. In the Appendices are included details of the derivation, 
which were omitted in the main text in order to facilitate the reading. 
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2. Theoretical Background 

We consider a system of localized effective spins characterized by the Hamiltonian 



^ = ^cxc + ^dip + + ^SR + ^-R + ^SL + ^L: 



where 



^exc — ^ ^ J(Rij)Si ' S 



ft: 



(2) 



accounts for the exchange interaction between pairs of localized spins Sj and Sj in 
the equilibrium positions Rj and Rj of the magnetic ions that are present. With 
more than one per unit cell indexes i and j would be composed of the one indicating 
the position of the unit cell and those indicating the positions of the magnetic ions 
within the given unitary cell relative to the cell position. The exchange integral J 



depends on the distance Rij 
interaction given by 



|Rj — Rj|. The second contribution is the dipolar 



dip 



E 



Si ■ Sj 3 (Si • Rjj)(Sj • Ry) 



R 3 - 



R 5 - 



(3) 



with g and /ie being the g-factor and Bohr's magneton respectively. The terms 
J^z, ^sr and J^r are related with the magnetic field present in the material: 



(4) 



is the so-called Zeeman term involving the coupling with an external constant 
magnetic field H ; the time-dependent magnetic fields H(t) (including the 
pumping rf-fields) are incorporated through the term 

J% R =-gV B tt(t)-J2Si> ( 5 ) 



while J#b. is the Hamiltonian of the free photons in the electromagnetic fields. 

The last terms account for the effects of lattice vibrations: is the 

Hamiltonian of the free phonons and the spin-lattice interaction is given by 

d 



SL 



E 



dK 



^s(Rij) + 



+ 



d 



(6) 
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with 9 rijt ^s(Rjj) evaluated on the equilibrium positions Ry and xy = Xj — Xj 
with Xj being the displacement of the ion around the equilibrium position Rj. 

Introducing a second quantization formalism (which, in the case of spins, is 
done in terms of Holstein-Primakov and Bogoliubov transformations [28| |29| 130] ) 
we arrive to the expression of the transformed Hamiltonian given in Appendix A 
[Eq. (1A.22j) ]. now in terms of magnon, phonon and photon creation (annihilation) 
operators: c q (c q ), 6 q (6 q ) and d q (d q ) respectively. We call the Hamiltonian of 
Eq. (lA.22j) the magnons' Hamiltonian, which we use in the calculation of the 
nonequilibrium magnon populations. 

For that purpose, first, the thermo-statistics deemed appropriate for the 
description of the nonequilibrium macroscopic state of the system needs be 
introduced. As noticed in the Introduction we resort to the use of NESEF. 
According to the formalism, following Mori, Zwanzig and others [3T| l32j 133} 13^1 
l35| 136] . the Hamiltonian of the system under consideration is separated out in a 
so-called relevant part ji?o, consisting of the energy operators of the free degrees of 
freedom, and containing the interactions among them and the coupling with 
external sources and reservoirs. In our case here 

= ^ ^qCqCq + ^ ^k&k&k + ^ ^P^A> ( 7 ) 
q k p 

that is, the Hamiltonian of the free magnons [of Eq. (lA.llj) ]. free phonons [only 
acoustic ones considered, Eq. f lA.15j) ] and free photons introduced by the sample 
black-body radiation [Eq. (IA.2 ![) ] . with quasi-momentum q, k and p, and energy 
fiWq, hQk and h( p . On the other hand, M" contains the contributions associated 
to the energies of interaction presented in Eqs. f)A.22[) . 

Next step in the application of the formalism consists in the choice of a basic 
set of variables that should characterize the macroscopic state of the system (the 
appropriate nonequilibrium thermodynamic state of the system |H| US])- At the 
microscopic level, for the phonons and photons are taken the Hamiltonians 
and and for the magnon system is introduced the free magnon Hamiltonian 
j£g and the magnetic moment density operator 



M(r)=^B^S 4 5(r-R 4 ). 
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After applying Holstein-Primakoff and Bogoliubov transformations, we do 
have that 

M x (r) = g^ — J^e^ (v^ + Wq ) (cl q + c q J , (9) 
q 

M,(r) = ig^ — ^ e iqr (v^ - Wq ) (cL q - c q ) , (10) 
q 

q,q' 

X { (w*W q / + W q W*,) C q C q / + U q U q / C_ q C q / + 

+KK' C q C -q' + Vtf W } > ( 1 1 ) 

where in M x and have been conserved only the linear contributions, iV is the 
number of sites and we recall that u n and v q are the coefficients in Bogoliubov 
transformation, given in Eq. ( 1A.13|) in Appendix A. 

Therefore, a priori, the set of basic microdynamical variables is composed of 

[^ 2 \ M(r), ^ L + t # R }, (12) 

that is, the free magnon Hamiltonian M'^' of Eq. ( lA.llj) . the magnetization M(r) 
of Eq. (JSJ) and + of Eq. (J7J). However, for the purposes of the present 
study, it is convenient to refine this basic set introducing each of the contributions 
in Jffa and in Eqs. [9] [10] and [11] namely 



°~q — C-qCqj', |o"q — C_ q C q j; 



{4 



c q ;^ L + ^R k (13) 



where Q ^ 0. 

It can be noticed that this set consists of components in reciprocal space 
of the single-magnon reduced density matrix (Wigner - von Neumann single- 
particle dynamical operator [16] 07] EES] ) composed of the diagonal elements 
A q = c q c q , the occupation number operators, and the non-diagonal ones Aq q = 
a q +Q/2^q-Q/ 2 w ith Q ^ 0. The latter, describing the local inhomogeneities of 
the occupations A^ q , are not relevant for the present problem once space-resolved 
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experiments are not considered at this point. It can be noticed that all the 
single-particle observables of the system are expressed in terms of the single- 
particle reduced density matrix [47j . Moreover, have been introduced the creation 
(<r q ) and annihilation (cr q ) operators of magnons pairs, and their nondiagonal 
contributions (associated to local inhomogeneities) are not considered for the same 
reason appointed above. Finally, being bosons, are introduced the creation and 
annihilation operators in magnon states, c q and c q , whose eigenstates are the 

~ (2) 

so-called coherent states. After this considerations and recalling that is 
expressed in terms of occupation number operators [see Eq. (lA.llj) ]. the set of 
microdynamical variables relevant to our problem is then composed of 

{{^} ; {^}; {*,}; {^}; -4 + (w) 

According to NESEF, the nonequilibrium statistical operator, given in 
Appendix B, depends on the quantities in set f|T4|) and on another set of 
nonequilibrium thermodynamic variables associated (also said thermo-dynamically 
conjugated) to the basic ones in (T14T) which we designate, respectively, by 

F^t)}; {*>;(*)}; {¥>„(*)}; {^(t)}; A> }, (15) 

with (3o = l/ksT [see Eq. ( IB.Gj) ]. Finally, the space of nonequilibrium 
thermodynamic variables consists of the average values over the nonequilibrium 
ensemble of the quantities in set (JHJ), say, 

|{jV q (t)}; |^ q (i)}; {(^l*)}? ^ + ^|, (16) 

that is, 

Af^t) = Ti^g e (t) x (17) 

with g £ (t) of Eq. (\B.2\i and of Eq. (IB.6p . and so on for all the microdynamical 
variables in the set of Eq. (JHJ). 

We are now in conditions to go over the derivation of the evolution equations 
for the set f lTB]) of basic macrovariables, i.e. to obtain the time evolution of the 
nonequilibrium thermodynamic state of the magnon system. 

3. Evolution of the Nonequilibrium Macrostate of the Magnon System 



The equations of evolution for these variables are the quantum mechanical 
equations of motion for the dynamical quantities of set ( Tl4|) averaged over the 
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nonequilibrium ensemble. They are handled resorting to the NESEF-based 
nonlinear quantum kinetic theory [31] HOI HTJ H2J S3], with the calculations 
performed in the approximation that incorporates only terms quadratic in the 
interaction strength, with memory and vertex renormalizations neglected |40[ W2\ 
I5T)] . that is, we keep what in kinetic theory is called the irreducible part of the 
two-particle collisions. In the case of the population of magnons the equations of 
evolution are given by 

d . , , s 1 ^ 



in 
+ 



jVq, 3tf 



Qe(t) X g B 

g(t, 0) x Qb\ + 



1 



[ihf 



dt' e'V-Q Tr 



M"it! -t)^'M^ g(t,0)xg B + 



5Q,(t) ' 



where g(t, 0) is an auxiliary statistical operator (cf. Appendix B), e — > 
after the calculation of averages, lower index naught indicates interaction 
representation, Q e stands for the quantities in the set of Eq. ([IE]) and jj^j = 



^)- lr Tr 



A/- q ,JT' 



g(t) >. In the case of amplitudes, 



d , . 1 m 

* < c < |f > = il 1 * 



-Tr 



Cq, 



Cq, 



Qe(t) X £ B 
g(t, 0) X ^b \ 4" 



+ [ dt> e£( *'^ )Tr { - *)o, c q ] Q(t, 0)xg B \ 



5Qi(t) 



(19) 



the one for (c q |t) is the complex conjugated of this Eq. ([19"]) . and the evolution 
equations for the magnon pairs are 



d , s l m 



+ 



Tr 



<7q, 



g(t, 0) X g B \ + 



[ihf 



dt' tfV-Q Tr 



- t) , [JT, <x q ] g(f , 0) x g B 
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+ ^ ee(t '^{[^ / - i )o'^]^0)x^ B }^|,(20) 

and its complex conjugate for cr*. 

They acquire the quite cumbersome expressions shown in Appendix C. Here 
we present them in a compact form, indicating and describing the contribution of 
the different processes involved, which for the populations is 

j t K(t) = © q (t) + *Kq(t) + L q (t) + £ q (t) + $ q (t) + 97t q (t) + 2t q (t) , (21) 

where on the right we do have: (i) © q (t) is the rate of growth of the population 
in q-mode produced by the external source, which is composed of 2 contributions, 
namely, a direct production, and a term of a positive feedback (only associated to 
parallel pumping excitation); (ii) £H q (£) is a non-linear term of relaxation due to the 
decay of magnon in photons, leading to the saturation of absorption when under 
continuous excitation; {Hi) L q (t) is a term of linear relaxation to the lattice with a 
relaxation time r q ; (iv) £ q (t) is a term involving nonlinear relaxation to the lattice, 
referred to as Livshits contribution|51J; (v) 3q(t) is a peculiar and fundamental 
contribution of a nonlinear character arising out of the magnon-lattice interaction 
[the sixth and seventh contributions in Eq. (1C.22I) of Appendix C], which takes 
the form 



2% 

qVq 



ft 2 



+ fi2 Yl {(-^ + ^-A/q'^q-q' " -A/qC-A/q' + l)(^q-q' + 1)} X 

x% q( -w q + n q _ q /), (22) 



where z/ q and f2 q are the population and the frequency dispersion relation of 
the phonons in the thermal bath [see Eq. (1C.22|) in Appendix C]. After some 
mathematical handling, this Eq. (122 p can be rewritten in the form 

3q(*) = E Xqq' {^q'(^q + 1) 6*^' - (A^ + l)jV q } , (23) 

q' 

where 

2tt 



Xqq' = |^q,q-q'| 2 {>q-q' e 5{u)qi - W q + Q q _ 



/l 2 



q' 



+ ^q'-qe-^^^K ~ " q - *Vq')}> (24) 
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with this Eq. ( 123|) having the form given originally by Frohlich |11[ IT2]. and we call 
it Frohlich contribution. 

Considering high population values {Nq> 3> 1), Eq. (1251) becomes 

Ut) = E Xrt^K (e A ^' - e*^), (25) 

q' 

and we can see that, since Xqq' > 0, the contributions for the Eq. fl25]) are positive 
for those modes q' for wich Uq> > a; q and negative for those modes with Uq> < u; q . 
Consequently, modes q' for which u)q> > w q transfer their energy in excess of 
equilibrium to the mode q, and therefore in a cascade-down process it is transferred 
to the mode lowest in frequency. Thus, the mode lowest in frequency largely grows 
in population (drained from all the other modes) leading to the emergence of, what 
has been dubbed, a nonequilibrium Bose-Einstein condensation. Moreover, we 
emphasize that the Frohlich term has a purely quantum mechanical origin. We can 
summarize the point stating that such nonequilibrium Bose-Einstein condensation 
of "hot magnons" is of a pure quantum character and driven by Frohlich non- 
linear contribution to the kinetic equations, whose origin is in the interaction with 
the thermal bath in which the system is embedded (a description of the irreversible 
thermodynamics involved is presented in Ref. [Hj). 

The other contributions in Eq. (1211) are: (vi) 9Jt q (t) is the rate of change 
generated by the magnon-magnon interaction (exchange and dipolar), whose role is 
to lead the system of magnons to a state of nonequilibrium internal thermalization. 
This contribution is in a "tug of war" with Frohlich contribution (previous 
item); (vii) 2l q (t) contains all the contributions coupling the populations to the 
amplitudes, (c q |t) and (c q |t), and the pair functions, o" q (t) and o"*(t). 

Therefore this evolution equation is coupled to those of these other basic 
variables given in Appendix C. The evolution equation for the amplitudes, also in 
a compact form, is given by 

j t (c q |t) = -zu; q (c q |t> - r q (t) (c q |t> , (26) 

where, in Mori's terminology [52], the first term on the right is a precession term and 
the second, cf. Eq. ( IC.19p . is in balance a relaxation (damping) term containing 
contributions arising out of the magnon-phonon interaction (of linear, Livshits 
and Frohlich type in the nomenclature already introduced), of interaction with 
the radiation fields, and from the magnon-magnon interaction. 

Similarly, for the pair magnon function follows [cf. Eq. ( 1C.20|) ] that 

-(7 q (t) = -2za; q a q - T q (t) a q + A q (t), (27) 
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which is the first harmonic of the basic one in Eq. (|26|) . and A q (£) consists of a 
term involving the effects of the external rf source and contributions coupling via 
interaction with the latter and the magnon-magnon interaction, with the other 
a q.'(t) Q.' 7^ Q- Solution of Eqs. (j2Tj) , f j26|) and f )27|) requires to provide initial 
conditions. Considering the initial state as a equilibrium one, the initial condition 
for populations A/q(t) is Af(t = 0) = (e^ ^ - l)~ , and those for (c q |t = 0) and 
cr q (£ = 0) are zero. Therefore, in the conditions to be analyzed, (c q |t) remains null 
(i.e., there is no contribution to the total population from the population of the 
coherent states). On the other hand, |cr q (t)| tends to increase with source A q (t) 
(and to decay with a lifetime given roughly by a time average of T~ 1 (t) during the 
length of the process), but the increase involves a rate of change similar to the one 
for jVq(t) and then, in quite general conditions, the contribution to the population 
due to the one of the magnon pairs is orders of magnitude smaller than the leading 
one corresponding to that of the individual quasi-particles (single magnons). 

This can be seen in the fact that, on the one hand, Eq. fT2B]) can be expressed 
in the integral form 

(c q |t) = (c q |0)exp|^ rft'H Wq (t')-r q (t')]}, (28) 
and, on the other hand, Eq. (T2T1) becomes 

<Tq(t) = <7 q (0) exp j jf dt' [-i2w q (t) - T^t')] j + 

+ exp j jf dt' [-i2w q (t) - r q (0]| x 

x jf dt' Aq(t') exp I jT dt" [i2w q (t) + r q (t")] I . (29) 

Moreover, a direct calculation provides us with the nonequilibrium 
thermodynamic equations of state, that is, the relation between the basic variables 
of the set (1161) with the nonequilibrium thermodynamic variables of set (fT5|) . 
namely 

A/- q (t) = Tr {A/q g £ (t)} = A/" q m (t) + A^ oh (t) + A^ air (t), (30) 
where we have introduced the definitions: 

^rw - ( 31 ) 

is the population of the single magnons; 



A/q C ° h (t) = |(c q |t>| 2 



0q(*) 



F^) + bq(t) + ^-q(*)] 



(32) 
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is the populations of the coherent states; and 

yv i W - e F q _ ! | 2F'(t) - 1 



e^-1 1 

is the population of the magnon pairs; F^(t) = yj Fj*(t) — |y q (i) + v?_ q (i)| 2 . See 
appendix D for details. It can be noticed that we can redefine the nonequilibrium 
thermodynamic variable -F q (£) in either of two ways, 

F q (t) = /3 [/^ q - /i q (t)] or F q (t) = huJk B T*(t) , (34) 

following, respectively, Frohlich [11 J and Landsberg [53J introducing a so- 
called quasi-chemical potential per mode, // q (£), and Landau, Uhlenbeck and 
others (a description in [3H IMj), introducing a so-called quasi-temperature 
(or nonequilibrium temperature) per mode [HH [55], T*(i), as it is usual in 
semiconductor physics EZ1 • 

As described above, A^ oh (t) = and A^ air (t) < A/^ m (t), and then in Eq. 
21] we can take 2l q (£) ~ (i.e., 2l q (t) orders of magnitude smaller than the other 
terms) and then Eq. [2T]is closed in itself and we proceed to deal with it. 



4. BEC in YIG 



For numerical calculations and comparison with experiment, first we introduce a 
treatment consisting in a kind of "two-fluid model", namely, we transform the 
large system of coupled evolution equations in a pair of coupled equations for the 
mean values of the populations over two regions of the Brillouin zone: one is a 
small region around the position of the minimum in frequency in the dispersion 
relation (~ 2, 1 GHz), which we call R±, and the other around the zone in which 
are the modes pumped by the rf-fields (~ 4 GHz), indicated by i? 2 - Such kind of 
procedure is justified, first, because of what followed in the other cases of BEC 
in bosons (phonons, excitons) where a complete solution is obtained (symmetry 
conditions allowed to separate the set of coupled equations in blocks with a small 
number of coupled equations) showing that kind of behaviour and, second, it is 
verified to a good degree in the experiments, as can be noticed in Fig. 3 of Ref. 
[1] and better in Fig. 2 of Ref. [2]. The chosen regions are shown on the right 
of Fig. [1] and were determined on the basis of what is shown by the experimental 
Brillouin spectra reported in Refs. [I] and [2]. 
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Figure 1. On left we see the thin film geometry, with magnons propagating in 
the rrz-plane. On right, the regions R\ and R2 in quasimomcntum space. 



4-1. Evolution of the "Two Fluid Model' 



In terms of the previous considerations we introduce the mean populations M.(£) 
and M 2 (t) corresponding to the relevant regions R\ and R 2 given by 



1 qei?i 



(35) 
(36) 



qei? 2 



where n i)2 = 1 represents the number of modes in the regions R\ and R2. 

qe/?i, 2 

Proceeding correspondingly in Eq. [2TJ we do have that (omitting to explicitly 
write the time dependence on the right) 



/i|-M(t) = -DM(M - A/f) - /1 [M - A/f 

+ F {MM + (u + l)Af 2 - pMi} - 

- M {M (M + 1) + M (M + 1)} (M 



+ 



A/f 
M (0) 



M), (37) 



and 



f 2 jM 2 (i) = l(l + 2N 2 ) 



dm 2 (M> - m 2 (0) ) - / 2 [m - M (0) ~ 

F {MM + (z/ + 1) M 2 - z>M} + 

M (0) 



+ m {M (M + 1) + M (M + 1)} (M 



M (0) 



M), (38) 



where t is the scaled time t/r, taking the relaxation time r q as having a unique 
constant value (q- independent), M 12 are the populations in equilibrium, and 
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fi and /2 the fractions of the Brillouin zone corresponding to the two regions 
in the two-fluid model. Moreover, the coefficients M and F are the coupling 
strengths associated to magnon-magnon interaction and to Frohlich contribution 
respectively, D is the one associated to decay with emission of photons, and v is 
an average population of the phonons; the Livshits term can be neglected. Finally, 
the parameter I is related to the rate of the rf-radiation field transferred to the spin 
system, whose absorption, as noticed, is reinforced by a positive feedback effect. 
All these coefficients are addimensional, being multiplied by the relaxation time r. 

Consider now the experiment in Ref. pQ. The fractions f\ and fi follows 
considering the size of regions R\ and R2, and, given the frequencies associated with 
these regions, the values of the mean populations in equilibrium are Af^ = 3 x 10 3 
and Af^ = 2 x 10 3 . An effective absorbed power of 10~ 2 W would correspond to 
I w 10~ 3 . On the other hand, an estimate of the M2 steady state value allows to 
evaluate that D « 10~ n . On the basis of a transient time for attaining internal 
thermalization near equilibrium of the order of 200 ns, it can be estimated that 
M < 10 -12 . Being r of the order of a few microseconds (and adopting r = 1 /is), 
we are left with F as the only open parameter. 

Using the experimental data [TJ, varying the values of the parameters around 
those given above and adjusting F, it follows the good agreement of theory and 
experiment shown in Fig. [2J 

We proceed to analyse the several processes leading to the increase of the 
populations in Fig. |2J that are described in Fig. |5J in terms of the rates of 
increase and decay, the contributions present on the right of equations ( 157)1 and 
( 155)1 . Multiplying these rates by the total number of magnon modes n~3x 10 14 , 
we obtain the rates related to the total number of magnons in regions R\ and i?2- 
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Figure 2. Evolution of the magnon population. Circles represent Demokritov's 
et al. data for the low energy magnon population [T], with the pumping 
being switched off after 1 /is. Solid lines show magnon populations of low and 
high energy, obtained after numerical integration of Eqs. ([57)1 and ([55} using 
the following parameters: f x = 3 X 1(T 6 ; f 2 = 3 x lO -4 ; 7V 1 (0) = 3 x 10 3 ; 
7V 2 (0) = 2 x 10 3 ; I = 8 x 1CT 4 ; D = 4 x 1CT 11 ; M = 3 x 1CT 14 ; F = 2 x 10~ 6 . 



Besides the rate of pumping from the external source (designated by I), are 
present the contribution arising out from Frohlich effect, which is a pumping term 
for Ai and a decay one for J\f 2 , and the contribution due to magnon-magnon 
interaction, which redistributes the pumped energy among the modes, tending 
to drive the system to nonequilibrium internal thermalization (non-scattering 
mechanisms, i.e., decay and emission are discarded for not being of relevance: 
energy and momentum conservation are impaired [28J). These two effects (whose 
composition is referred as F+M) are responsible for internal interaction of magnons 
Afi and Af 2 . Two other contributions correspond to decay, namely, a linear term 
of decay to the lattice (L) and the bilinear one of decay by photon emission (D). 
The sign (+) for pumping and (— ) for decay is indicated in the inset. 

During the time interval when the pumping is on, the source creates magnons 
on region R 2 , while internal interactions (Frohlich and magnon-magnon) annihilate 
them. These internal interactions are responsible for creating magnons in region 
Ri, that decay mainly through photon emission. Thus, the nonequilibrium two- 
fluid system has the energy pumped on region R 2 , transferred to region R\ via the 
composition of Frohlich and magnon-magnon effects and while being lost through 
photon emission. 

The interplay of all these processes changes as the rate of the pumping source 
is changed. The results reported above are for the scaled rate I = 8 x 10~ 4 /is -1 . 
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Figure 3. (Color online) Contributions (scaled) from all the distinct processes 
to the evolution of N\ and A/a- With the same parameters used in Fig. [2j we 
plotted the rates of decay (— ) and increase (+) of A/2 (top) and A/i (bottom). 



4-2. The Steady State 

Using the same parameters but considering constant application of the pumping 
source, the steady-state populations are obtained as a function of the source scaled 
rate of pumping, what is shown in Fig. HI It can be noticed the existence 
of a critical pumping scaled rate (better saying, a rate threshold) after which 
there follows a steep increase in the population of the mode lowest in frequency, 
characterized by A/i, corresponding to the emergence of BEC. With increasing 
pumping intensity a second critical rate (rate threshold) is evidenced such that for 
higher values of I is observed internal thermalization of the magnons which acquire 
a common quasi-temperature [cf. Eq. (|34[) ]. This implies that the magnon-magnon 
interaction overcomes Frohlich contribution and BEC is impaired. 
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Figure 4. Steady-state magnon populations [solutions of Eqs. (|57l) and (I55l) ] 
as a function of the scaled rate of pumping due to the source, using the same 
parameters as in Figure [5J It can be noticed the existence of a window for the 
emergence of BEC, which follows at a certain threshold of intensity. In the linear 
regime, at low intensity, no particular complex behaviour follows, and at high 
levels of intensity magnon-magnon interaction overcomes Frohlich contribution 
and there follows internal thermalization. 
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Figure 5. Steady-state magnon populations as a function of the scaled rate of 
pumping when Frohlich term is disregarded (F = 0). 
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Figure 6. Steady-state magnon populations as a function of the scaled rate of 
pumping when magnon- magnon interaction is disregarded (M = 0). 

We can better appreciate the role of both types of interactions in figures |5] 
and IHl In Fig. the Frohlich contribution is "switched off" (F = 0), and the 
action of the magnon-magnon leading to internal thermalization (for I > 3 x 10~ 4 ) 
is evidenced. 

On the contrary, in Fig. O where the magnon-magnon interaction is "switched 
off" (M = 0), the emergence of NEFBEC (for I > 10 -4 ) does follow unimpeded by 
the magnon-magnon interaction. 

Returning to Fig. HJ let us consider the interplay of the several mechanisms 
that lead to the formation of the steady state, when the rates of change are 
balanced. The rate of change associated with all the mentioned processes, in 
the steady state, for a range of values of the (scaled) rate of pumping from the 
external source, are shown in figure in the upper part for J\f 2 and in the lower 
part for A/i, using the same notation of Fig. [3j It can be noticed the presence 
of three regimes in correspondence with those indicated in Fig. |H In all cases 
energy is fed to the system through the "A2 magnons" pumped by the source, 
while Frohlich and magnon-magnon redistribute it, creating "A/i magnons". In 
the range of scaled pumping rates up to roughly 6 x 10~ 5 , the linear regime, both 
A/i and A/2 magnons are relaxing predominantly through the linear decay to the 
lattice. For greater intensities the linear decay to the lattice is not sufficient to 
wholly absorb the system energy in order to maintain the steady state. Thus, it 
is observed an increase on the A/i populations, characterizing the condensate, and 
in order to maintain the system stationary other relaxation mechanisms become 
predominant: "A/i magnons" decay through photon emission and "A2 magnons" 
through the F + M term. Finally, for I > 2 x 10 -3 , the non-linear photon emission 
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decay gains relevance, the J\[% populations increase, and the internal thermalization 
regime is attained. 
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Figure 7. (Color online) Role of contributions (scaled) from all the distinct 
processes to the stationary state of Af\ and A2 . The rates for A2 (top) and Ai 
(bottom), with the same notation of Figure [3j are presented. 

All this interplay explains, in a sense, the transition between the three 
regimes: In the linear regime, as the rate of energy transfer from the source 
increases, we have that, at low scaled rates of pumping, the populations attain 
values near the equilibrium ones, the usual linear relaxation to the lattice, with a 
characteristic relaxation time, predominates and no particular complex behavior 
is to be expected. As the intensity approaches a threshold value, the nonlinear 
contributions, associated to Frohlich effect, magnon-magnon interactions and 
photon emission, begin to be relevant as the populations increase to values much 
greater than in equilibrium. Frohlich effect leads to the emergence of BEC once 
it overcomes the tendency for internal thermalization promoted by the magnon- 
magnon interaction, and the decay by photon emission ensures that a steady state 
is attained. The second threshold in intensity value follows as magnon-magnon 
interaction overcomes Frohlich effect and internal thermalization is ensured. 

Moreover, considering the quantity T q of Eqs. [26] and [23, in this two-fluid 
model we do have for Ti and T 2 (mean values of T q in regions R\ and R2) that 

ri(f) = D (M - A/f) + ft - F (jV 2 - P) + 
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+ M 



A/f 



[M (M + l) + M 2 (A/2 + i)] - (M + 1) A/2 



21 

(39) 



and 

r a (f) 



21 + D (A/ 2 - A/" 2 (0) ) + / 2 + F (A/1 + z> + 1) + 
- I\ I <! AT (M + 1) + A/" 2 (A/" 2 + 1) - -^yM (A/ 2 + 1) 



(40) 



which have their stationary values shown in Fig. [BJ as a function of the rate of 
pumping I. 




Scaled Rate of Pumping I 

Figure 8. Values of I\2 in the steady state (cf. Eqs. l5§leH0]) associated with 
the magnon populations (Figure 0]). 



According to the considerations presented in section [3j we see that T q is the 
linear decay rate of the amplitudes, (c q |t) and (c q |t). All the obtained values of 
ri j2 are positive, (Fig. [8]), indicating decay of these amplitudes, which corroborates 
the assumed neglect of their contribution to the populations' evolution. Although 
the amplitudes decay, we point that in the source intensity interval associated with 
the condensate, 1CT 4 <I<2x 1CT 3 , Y x is considerably smaller than r 2 , and the 
condensate would allow long mean-life to the coherent states corresponding to the 
lowest frequency modes, which can be excited to compose solitary waves, as it has 
been shown in the case of other boson systems displaying NEFBEC [T4"l [TSl 20J . 

The decay of the populations after the switching off of the pumping source 
follows in accord with the one observed in the experiment of Ref. [SB]. It consists 
of three regimes: a near exponential one at the initial delay times (scaled time in 
the interval 1 to ~ 1.7 in Fig. [2]), with scaled decay time of 0.23, and another with 
a scaled decay time of 1.02 when approaching final equilibrium (after scaled time 
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~ 2.5 in Fig. [2]) and an intermediate one in between (roughly the interval from 
1.7 to 2.5 in scaled time in Fig. [2]) as shown in next subsection. 

4-3. Decay of the Condensate 

In Ref. [58] is reported by Demidov et al. an analysis of the decay, towards final 
equilibrium, of the NEBEC in YIG, after the external rf-pumping source has been 
switched off. There are some differences in the experimental protocol with respect 
to the one used in Ref. [1] , but we analyse the decay in the conditions of the latter, 
to study the dependence (as done in the experiment of Ref. [SB]) with the source 
power. 

Using Eqs. f!57|) and in Subsection 14.14 the one that leads to the results 
of Fig. [2J varying the values of the rate of pumping I we obtain the set of curves 
for the evolution of the population in the condensate shown in Fig. [9j 




Time (/is) 



Figure 9. (Color online) Time evolution of the condensate (Afi) for different 
values of the rate of pumping I. 

It may be noticed the interesting result that these curves can be approximately 
well adjusted by the law 

- M 0) = A exp (-*/ TA ) + B exp (-*/ T ) , (41) 

The fitting values of the two coefficients, A and B, and the characteristic 
decay time ta are indicated in Table 1; r is the relaxation time to the lattice. 
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Table 1. Coefficients A, B and ta obtained from the fitting of Eq. @J]to the 
numerical data presented in Fig. [5] (it has been taken r = 1 /j,s). 



I 


ta (a*s) 


A 


B 


5 x 10" 5 


0,47670 


140, 10 


1064,48 


6 x 10" 5 


0,46853 


217,96 


1316,59 


8 x 10" 5 


0,45151 


458,56 


1869,11 


10 x 10" 5 


0,4326 


863,2 


2496, 7 


13 x 10" 5 


0,3990 


2006,9 


3618,7 


16 x 10" 5 


0,3568 


4460 


5048 


20 x 10" 5 


0, 2849 


1,3603 x 10 4 


7786 


25 x 10" 5 


0,2131 


5,890 x 10 4 


1,278 x 10 4 


32 x 10" 5 


0,21706 


1,5904 x 10 5 


1,523 x 10 4 


40 x 10" 5 


0,22634 


2,3599 x 10 5 


1,524 x 10 4 


50 x 10" 5 


0,23113 


3,0488 x 10 5 


1,517 x 10 4 


63 x 10" 5 


0,23310 


3, 7937 x 10 5 


1,544 x 10 4 


79 x 10" 5 


0,23263 


4, 6423 x 10 5 


1,636 x 10 4 


100 x lO" 5 


0, 22948 


5, 765 x 10 5 


1,852 x 10 4 


126 x 10~ 5 


0,22314 


7,315 x 10 5 


2,29 x 10 4 


158 x 10~ 5 


0,2124 


9, 765 x 10 5 


3, 17 x 10 4 



Clearly, we can say that there exist three reasonably well defined regimes: 
one immediately after the switching-off of the source, other at longer decay times 
when the population is approaching the final value at equilibrium, and, of course, 
an intermediate one. 

In Fig. [10] we indicate the decay time ta as a function of the scaled rate of 
pumping I. 

It can be noticed a qualitative agreement with the experimental data of 
Demidov et al. in figure 3 of reference [58J, namely, a sigmoid-like curve up 
to I ~ 1 x 10~ 3 , except for the central dip in an intermediate region where 
no experimental data are reported. We proceed to study in detail the different 
relaxation processes that produce these results. 

We analyse (using the same notation for rates of the previous section) the 
relaxation processes involved: (i) the linear relaxation to the lattice (L-process); 
(ii) the radiative decay (D-process); (in) the conjugated effect of magnon-magnon 
interaction and Frohlich effect [(F + M)-process]. 

At small I, as stated before, the L-process dominates the decay. But, for 
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Figure 10. Decay time ta of the condensate as a function of the pumping 
parameter I. 

intensities such that Eq. HH applies to the populations decay (see figure [H]), it 
can be noticed the increasing relevance of the F + M-process, as can be seen in 
Fig. [TTr (Regime I). While the pumping source is acting, Frohlich and magnon- 
magnon terms transfer the energy from the fed magnons (A/2) to the ones in the 
condensate (A/i), but after turning off the source the flux of energy is inverted, and 
F + M-processes become another mechanism of relaxation of the condensate, thus 
reducing ta- Progressive amplification of the pumping power enlarge this effect, as 
shown in Fig. [TU] A minimum ta value is achieved for I ~ 2, 5 x 10~ 4 (Fig. [TTb) 
and then, for higher I, the F + M-process, dominant relaxation process until that 
intensity, diminish, and the D-process begins to increase. Then, the decay time ta 
increases until the D-process dominates the relaxation (Fig. [TTb ) and, after this 
point, where is initiated the Regime III, ta starts to decrease again. For values of 
I higher than 2 x 10~ 3 , the relaxation in the condensate does not follow Eq. HTJ 
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Figure 11. (Color online) Scaled rates of creation and annihilation of magnons 
in R\ region for different scaled pumping rates. On the upper figure, with 
I = 10~ 4 , in Regime I, the linear relaxation is predominant, but the F + M- 
process gains relevance. On the central figure, in Regime II, with I = 2, 5 x 10~ 4 , 
and the F + M-process is the principal relaxation mechanism. Finally, on 
the lower figure, we have in Regime III, with I = 10~ 3 , where the D -process 
dominates the relaxation. 



Thus we conclude that, depending on the populations A/i and A/2, immediately 
after turning off the pumping source, the three distinct contributions indicated rule 
the condensate relaxation. Briefly, it can be stated that increasing the pumping 
source intensity before turning it off changes the dominant relaxation process in 
the following order: 

L^F + M^D, 

explaining in this way the unexpected behaviour of A/i while returning to 
equilibrium. 
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5. Concluding Remarks 

In summary, as experimentally evidenced by Demokritov et al. the spin system 
in magnetic thin films under excitation by rf-radiation, display a phenomenon of 
the type of a Bose-Einstein condensation of the hot magnons pQ. The theoretical 
analysis performed here, in terms of a nonequilibrium statistical thermo-mechanics, 
shows that, in fact, this BEC of magnons follows in the way expected for many- 
boson systems embedded in a thermal bath as firstly demonstrated by H. Frohlich 
[HJ [12], and we may cite him who, in his original work (the case of biopolymers), 
stated that 

"... under appropriate conditions a phenomenon quite similar to Bose 
condensation may occur in substances which possess longitudinal electric 
modes. If energy is fed into these modes and thence transferred to other 
degrees of freedom of the substance then a stationary state will be reached 
in which the energy content of the electric modes is larger than in thermal 
equilibrium. This excess energy is found to be channelled into a single 
mode - exactly as in Bose condensation - provided the energy supply 
exceeds a critical value. Under these circumstances a random supply of 
energy is thus not completely thermalized but partly used in maintaining 
a coherent electric wave in the substance." 

Moreover, that long-lived solitons (in the form of a statistically averaged 
macroscopic wave function of coherent states) propagate in the condensate, Finally, 
it can be noticed that the extended analysis of the nonequilibrium thermodynamics 
of the phenomenon described in the Refs. [H] is formally identical to the one that 
can be applied to the magnons embedded in the lattice we have considered here, 
which is to be extensively described in a forthcoming article. 
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Appendix A. The Hamiltonian of the Magnon System 

To introduce the magnon creation (annihilation) operators, we use the Holstein- 
Primakoff transformation [28j [29j [30] : at first, spin operators are expressed in 
terms of local bosonic creation (annihilation) operators, dj (tZj), 



2 
2S 



2i 



a^aj 
~2S~ 



djdj 
~2S~ 



1/2 



a, + d]\ 1 



1/2 



a, — a] I 1 



2S 



1/2- 



S j — S dj CL j . 



(a.i; 



We consider a crystaline material with N c unit cells, with the ionic position Rj 
given by the position r n of the unit cell n and the internal relative position cLj, 

R = r n + d M , (A.2) 

with \l designating the internal ions (// = 1,2, .. . until the number of magnetic 
ions included in the unit cell), and it is introduced the Fourier expansion, 



e i(q-Rj) 



q,A»> 



(A.3) 



being q the wave vector running over the first Brillouin zone. Collecting only the 



quadratic terms of + J#di P + [Eqs. 
one obtains 



and (HD1 in what we call J$? G 



(2) 



S ' 



+ 2 a q>M a -qX 



(A.4) 



with 



A ^'(q) = E { Ai ( r -' + d^)A^' +d ^h 

n,n' 

+A 2 (rw + d^)} , 

B ^(q) = E B ( r -' + d ^)e i [ q(r -'- 



/+d„ 



(A.5) 
(A.6) 
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where r nn > = r n r — r n and with (n, //) 7^ (n f ,f/), and 

Ai(Ry) = -25J(i2y) - ^P[24 - 3 (^) 2 ]' ( A - 7 ) 

A a (R^) = 2SJ(J^) + ^ - ^^[4 - 3(i$) a ], (A.8) 

B(Rij) = _!(«fi)!scsj!_ (A 9) 

The diagonalization of J£g of Eq. IA.4l is done through the introduction of the 
magnon creation and annihilation operators c q7 and c q/r , as linear combinations 
of a„ and a q Ai in such a way that [2SJ ESJ ES] 

^ s (2) = ^ fru; q , 7 c qi7 c q , 7 . (A.10) 
q>7 

The resulting magnons, with energy ftw^ and group velocity V q u; qi7 , are 
grouped in branches indicated by 7, in the case of more than one magnetic ion 
per unitary cell (see Refs. [651 166] for the YIG). In the experiments with thin 
films of YIG analyzed here, only low frequency magnons are excited (< 10 GHz), 
justifying the omission of all the branches but the acoustic one. In this sense, only 
one effective spin per unit cell is considered, fi and 7 are omitted and we write for 
J^ 2) of Eq. flAlOD 

Jg£ 2) ~ ^qCqCq- (A. 11) 

q 

In this situation of taking into account only the acoustic magnons, Eq. (lA.lip 
follows from Eq. (IA.4I) . when in the latter we take one effective spin per unit cell, 
after using the so-called Bogoliubov transformation, [28| 129] , 

* t 

fl q W q C q T" V „C_q, 

fl -q = M -q C -q + W -q C q> ( A - 12 ) 

with Uq and v„ functions of q. It can be shown that 



/ A'(q) + ^ q _ B'(q) A(q)-^ q 

Mq -\/ 2^ q ' Uq -pMV ' ( } 

and the (acoustic) magnon dispersion relation is given by 



^q= V[ A '(q)] 2 -|B'(q)| 2 , (A.14) 

and Refs. (H |67] presents recent studies on this dispersion relation in cases of thin 
films of YIG. 
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The Holstein-Primakoff and Bogoliubov transformations are then applied to 
the non-quadratic terms of J££ X c + <^<iip + ^z, and the magnon-magnon interaction 
term is then obtained. Retaining only the fourth order scattering terms, the 
magnon-magnon interaction which contributes to the Hamiltonian is given by 

^MM = ^ ] 1 /? q,qi,q2^q^qi^q2^q+qi-q2- 

q>qi,q2 

Phonon and photon creation and annihilation operators are introduced in 
similar ways. The Hamiltonian of the free phonons is 

j#L = ^m k (htb k + ~^ , (A.i5) 

with f2k being their dispersion relation (we recall that only acoustic phonons 
were considered; polarization is implicit). The phonon creation and annihilation 
operators, b\, and 6k, are related to the displacement of the effective magnetic ion 
around the equilibrium position r n , which is 

--(^)'?^(^ + ^"' k "')- (A16) 

and e(k) is the polarization versor (cf. Ref. [2H])- Using Eq. IA. 161 in Eq. El we 
obtain the for magnon-phonon interaction 



q,k^0 



^sl = (&k + 6^ k ){Vq,kC q c q _ k + £ qik c q 4_ q + 

+ £q,- k CqC-k-q} + 

q,k^0 

+ ^-;,-k^k6k-q}(Cq + C t _ q ), (A.17) 



where .Fq,k 5 £ q ,k an d ^• q ± k are the resulting coefficients (representing the 
interaction coupling intensities). We write for the field generated by photons of 
the electromagnetic fields (from the source and black-body radiation) [68] 

H(r) = ^(H QiP (r) d a , p + H* iR (r) ^ p ), (A.18) 

ct,p 

where J aiq (d\ ) are the photon creation (annihilation) operators, 



H Q , p (r) = zp x A a>p (r), A Q , p (r) = J^ e ipr e^, (A.19) 
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C p is the photon angular frequency, p its linear moment and the polarization 
vector (the different polarizations are indexed by a). 

After some algebra, we obtain for the magnon-photon interaction of Eq. Othe 
expression 

^sr = £)(4p + dl_ p ) (S^cl + <Si_ p c_ p ) + 



q-p + ^q^q-p + 



+ ^a b q,pCqC P _ q + iS>jj, _ q _ p C_ q C q _ p [, (A. 20) 



«.p>q 

' '•i.q-p' q' p-q <-M.-q.-p' -q' q-p J- 

where 5^ P , Sa%, P and iS>l b q )P are intensity coupling coefficients. Considering that 
the external pumping source and black-body radiation are the origins of the 
electromagnetic fields, we can distinguish respectively their operators as d| q and 
<\^ q , and the energy of the black-body radiation is 



^ = E *Cp feiC + I) ■ (A.21) 

o,p ^ ' 

The complete Hamiltonian is then written, in second quantization form, as 

J? = ^ ^qCqCq + ^ ^k&k&k + ^ ^Cp^,q^a,q + 
q k a,p 

+ ^ ] "' / q,qi,q2CqCq 1 Cq 2 Cq+q 1 _q 2 + 

q,qi,q2 

+ J2(b k + 6 f _ k ) {-Fq.kCqCq-k + >Cq,kCq4-q + ^-k^qC-k-q } + 
q,k^0 

+ J2 {^q^k- q + 7^ k ^ 

q,k^0 

+ J2(<p + 4- P ) ( 5 ipcp + <_ p C- P ) + 

a,p 

+ 5^ (^ a 'P ^a,-p) |^I,q,p^q^q-P + ^l,q,p^q^p-q + ^Q,q,-p^-q^q-p | ■ 

",p,q 

(A.22) 

Appendix B. The Nonequilibrium Statistical Operator 

The nonequilibrium statistical operator is given by [3TJ [32j [33j EH ESI EE] 

= Qe(t) X §B. (B- 1 ) 
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where 

Q s (t) = exp |ln^,0) - ^W^^ln^t' - *) j (B.2) 

is the nonequilibrium statistical operator of the magnon system, with the auxiliary 
statistical operator p (also called "instantaneous quasi-equilibrium operator"), 
depending (superoperator) on the basic microvariables of set [HJ given by 

p(t,0)= expj-$(t)-^[F q (t)c q c q + 
^ q 
+ q (£)c q + q (£)c q + 

+^ q (t)c q c_ q + ^(t)c q cl q ]| (B.3) 

after recalling that we haze neglected local inhomogeneities. In p(t', t")the first 
term in the argument, t', refers to the evolution in time of the nonequilibrium 
thermodynamic state of the system (i.e., of the nonequilibrium thermodynamic 
variables Fq(t), q (t) and </? q (t)), and the second, t", to the evolution of the 
microdynamical variables A/" q , c q and <r q in Heisenberg representation, e is a 
positive infinitesimal that goes to +0 after calculation of average values has 
been performed. Moreover, ensures the normalization of the probability 

distributions, and plays the role of the logarithm of a nonequilibrium partition 
function: $(£) = lnZ(t). It is verified that 

*•,<«) - |g, (B.4) 

and similarly for the other quantities, in complete analogy with the situation in 
equilibrium, meaning that -F q (£) and the others in set ( !T5|) are the nonequilibrium 
thermodynamic variables said conjugated to the basic variables in set ( II 6p . 
Furthermore, it is verified that 

g £ (t)=p(t,0) + g' e (t), (B.5) 

where g' e (t) incorporates irreversibility and historicity. In Eq. IB. II 

& = ^ B 1 (To)exp|^±pJ (B.6) 

is the canonical distribution function of the phonons and photons in equilibrium 
at temperature T . 
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Appendix C. The Evolution Equations for the Basic Variables 

The macrovariable Qj [a generic expression for those of the set indicated in Eq. 
(flBI) ]. related to the microdynamical one Pj [of the set in Eq. (fill) ], has its evolution 
given by the Heisenberg equation of motion weighted with the nonequilibrium 
statistical operator M £ {t), namely 



s*« = if 1 



d 



with 



4>m 



1 
1 



Tr 
Tr 



Pj, JT° 



g{t,0) x g B 



Pj,J#" g(t,0)xg B 



and 



Pj,J?" 



g' s (t,0) x £ B } 
In the Markovian approximation we have that 

dt'e^'-^x 



(my 



x Till 3fc\H - t) , Pj] g-(t, 0)xg B \ + 



+4z 



dt' e £(i '^ x 



5Qe(t) 



x Tr \[j?'(t' -t) ,Pj) g-(t,0) x g B 



(CI) 



which, considering Eq. (IB.lj) and that, according to \B.5\ g' £ {t) = g e {t) — g(t,0), 
can be rewritten as 



(C.2) 

(C.3) 
(C.4) 

(C.5) 



(C.6) 



and the evolution of Qj, is thus expressed only in terms of average values weighted 
with g(t,0), where S stands for functional derivative. 
The evolution equations for the amplitudes are 

° ' ' " ■ (C.7) 



dt 



(^\t) = j£\t) + j£\t) + j£\t), 
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with 

Jg ) (*) = -^q<«ql*> 

(the precession term in Mori's terminology [52J), 



},qi,q 2 \^qi C q2 C q+q 



i-q2 1^) 



qi,q2 



z /i 



q+qi -q 2 



qi,q2 



qi 



qi 



arising out of the magnon-magnon interaction, and 
J^t) = J?){t) l + J?)(t) ll , 

where 



^ Q 2) (t)i = J { 2(t)f M + 4? (Of + 4l\t)! 



SR. 



with 



and 



x Tr{ \je"(t'-t) ,Pe Q(t,0) x g B 



x Tr 



^VUVl(£' - Oo, [^VIM, C q ] £(t, 0) X £ B j 



x Tr 



(my 



^Sh(t' - t) , [J&SL, Cq] £(t, 0) X Q B 



dt' e<*'-*) x 



x Tr 



^SR^' - t) Q , [J%R, Cq] g(t, 0) X ^ B } • 
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(2) (2) 

The scattering integrals Jc q (t)i and J Cq (£)n should be written in terms of 
populations, amplitudes and pairs of magnons, and in which double commutator 
of Eqs. (IC.12j) - (IC.14j) and the i-sum of Eq. fIC.llj) generate a huge number of 



terms, of which we analyze a particular one for illustration, say 

JcJ{t)l = — 8h ^ ] 1 /? q,qi,q2"' / 'q3,q2,q4 (^q 3 ^qi^q4^q3+q2-q4^q+qi-q2 1^) 

qi,q2,q3,q4 
J — oo 

+ Ah ^ ^ 1 /? q,qi,q2^q3,q4,qi (^q3^c|4^q2^q3+q4-qi^q+qi-q2 1^) 



X 



qi,q2,q3,q4 




x 



4^ ^ ] ^q,qi,q2^q2,q+qi-q2,q3 (^qi^q3^q+qi-q3 1^) x 
,q2,q3 

I (i 1 - e [ £+i K+ u q+qrq2- u i3- u q+ir'o)] T ) (C.15) 
J — oo 



qi,q2,q3 

X 



where r = t' — t. Average values are calculated and then expressed in terms 
of populations, amplitudes and pairs of magnons, as discussed in the following 
Appendix D. The average value enclosed in the last line, for example, gives 

(^qi^q3 Cq+ qi -q3 

I*) = [^(t) - |(6 qi |t)| 2 ] (5 q3 , qi (Cqlt) +5 

q3.q \ c qi 

+ [Vqait) - (Cq 3 \t) (C-q 3 \t)} 5 qi , q (c\\t) + 

+ (Cl\t) (Cq 3 |t) (Cq +qi -q 3 |t) . (G.16) 

The integration in time together with the limit of e — > +0 produces the so- 
called retarded Heisenberg delta function, that is 

1 



lim / dre [£+tullT = PV- -iir6(u), (C.17) 
where PV stands for principal value. 

(2) 

Applying this procedure to the all parts of J Cq one obtains the evolution 
equation of the amplitudes, Eq. (10. 7[) . Despite the extension of the obtained 
expressions, it is easy to note that all terms possess linear dependence on the 
amplitudes, as shown in the Eqs. (IC.8j) and (1C.8|) and because of the odd number 



of creation or annihilation operators in the commutators of Eqs. flC.lip - flC.14p . 
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Considering then a linear approximation and neglecting the self energy correction 
to the frequencies o; q we may write 

J t (Cq|*) = -l'Wq (Cq\t) ~ T q (t) (c q |t) , 

with frequency Co> q and 

r q (i) = Svr/t- 2 |V q , qi , q2 | 2 (7V q2 + AA q+qi _ q2 + l)AA qi x 



(C.18) 



qi,q2,q3 



-8?r/l 2 ^ l |/ q,qi,q2r- /V q2- /V q+qi-q2 X 

qi,q2,q3 

X <5(w q + W qi - W q2 - W q+qi _ q2 ) + 

+ 7T^ 2 ^ l-^q.kl 2 (A/q-k + + 1) <5(^k + ^q-k ~ Wq) 
k^O 



X 5(ujq +W qi -W q2 -( 

I"^q,qi,q2l A/q 2 A/q+ qj - -qj 



(A/q_k - Z'-k) <K^-k - ^q-k + W q ) 



+ 



+ 4vrn 2 ^ l^q,k| 2 (A/k-q - fk) <K^k - W q - w k-q) + 



k^O 



+ n 2 ^ |^-q,k| 2 (^k - ^k+q) £(^k+q ~ ^k ~ W q ) + 



k^O 



I 

+ 2fT 2 ^ ^ q , k (^-k + ^q+k + 1) <^-k + ^q+k " Wq) + 



k^O 

+ 7Tfl- 2 |5fqp| 2 (AA q _ p + / p + l) ( 5(Cp + 

a,p^0 

-Tt/T 2 ^ <5^q,p (JVq-p - /_p)5(C-p - W q _p + W q ) + 

+ 4irh- 2 |5ii b q, P | 2 (AA p _q-/p)<5(Cp-a;q-^p_q). (C.19) 

a,p^0 

The pairs and populations of magnons equations of evolution are obtained 
in an analogous proceeding. After calculating the collision integrals the evolution 
equation of the pairs of magnons may be written as 



^Cr q (t) = -2iw q Cr q - r q (t) CT q + Aq(t). 



(C.20) 



The first therm is associated with the frequency of precession of pairs, 2a; q , and 
the second is a decay term ruled by the same function r q (£) (Eq. IC.19I) related 
with the amplitude decay. A q (i) represent the non-linear terms, 

+ A^(t) + B^{t) (C.21) 



qW dt 



SL 



dt 
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da q 

The terms 



dt 



SL 



dt 



are originated in spin-lattice and spin-radiation 



SK 



interactions and depend exclusively on magnons populations; (t) are non-linear 
combinations of pairs and populations; the last therm represent the amplitude 
contributions to the pairs evolution. 

Finally, the populations evolution equation is 



dt 



8tt 



E K" q+q ^{(i + ^ q + AA q o/ q s , +q }5K + ^-C q+q O + 



h 2 



[©,(*)] 



qV-q 



+ ni E K 



b 

q>q + q ' 



q Y-q 



±[*,-jw 



(0) 

q 



+ 



{(AA q , + lXA/q + l)/J +q - A/q'A/"q(/ q / +q + 1)} Sfa + Wq , ~ C q+q " [<Kq(*)] 



[iq(*)] 



8tt 



+ -^2 E l £ qq+q'l {C^tf + + ^q'+q - A/q'A/q^q'+q + 1)} 5(w q + 0; q / - n q+ q') + [£q(t)] 



qV-q 



2tt 



+ TJ E l^qq-q'l {^'(A/q + l)(V-q + X ) ~ + iWq^'-q} *(Wq/ "Wq- O q ,_ q ) + 



h 2 

2ir 



q '^ q 



+-^2 E l-^qq-q'l {C^i + ^A/q'^q-q' - A/ q (7V q / + l)(i/ q _ q / + l)} $(c«v - ^ q + O q _ q /)+ 



[&,(*)] 



{(Aq + l)(AA qi + l)A/q 2 A/q 3 - A/qA/ qi (A/ q2 + 1)(A^ + 1)} X 



+ V" l v q>qi,q 2 l 

^ qi ^ q3 x5 K + w q i - - ^ q3 ) 5 q3 , q + q i- q 2 



+ E c i' ((^q'l*) ' (^q'lv ' CJ q '' Cr q '' t ) ' 



[2tq(t)] 

(C.22) 



where / p s and / p t are the average photon populations due to the source and 
black-body radiation respectively, and C q / ^(c q '|i) , ^c q ,|t^ , er q /, cr*,, are the 
contributions from amplitudes and pairs to the evolution. In compact form, as 
inEq. (ED, 

d 



dt 



A/q(t) = 6 q (t)+01 q (t)+L q (t)+£ q (t)+^ q (t)+9H q (t)+2l q (t) .(C.23) 
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The macrovariables of set in Eq. ( Tl6[) are the average of the set of microvariables 
in Eq. (Tl4l) weighted with the auxiliary statistical operator 



Q(t,0) 



exp I [^q^CqCq + 0q(*)Cq + 0* + ^ q (t)c q C_ q + (^ q (t)c q C f _ c 

I q 

Tr exp i -J2 [^ q (i)c q c q + + + ^ q (*)c q c_ q + </? q (t)c q ct c 



(D.l) 



This averages can be done through diagonalization of p(t, 0). Using the following 
transformation 



C q — t q ^ q + Kq^2 q + A q , 



(D.2) 



°— q q 1 "-q ^q i "— q) 

where and ^ q are bosonic operators and A q , 6 q and K q functions of q such that 

1 = [Cq,C q ] = 



^ q ^ q + ft q ^-cp ^q^q + K*^_ q 



^q, ^q 



+ kc 



- q) ^-q 



v q| 5 



it can be shown that 



(D.3) 



g(t,Q) 



exp I 

I q 

Trexpj-^F q (t)# q V q 



(D.4) 



0q(*) 



^q(*)+[^(*)+V-q(* 



(D.5) 



and 



1 F^t) + F&) 



2F>(t) 



(D.6) 



(D.7) 
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Using this transformation we show that the amplitudes are 
(c q |t) = Tr {c q §(t, 0)} = Tr | (i^ + Kq *fl q + A q ) g(t, 0)} = 

= ^ q Tr { <*f q 0) } + Kq Tr |Vl q 0) } + A q Tr { §(t, 0) } = A q , (D.8) 

(c q |t) =Tr{c q ^,0)}=Tr{(^ q t + K q ^L q + A q ) £(*,())} = 

= 6 q Tr §(t, 0)} + Kq Tr {^_ q |(f, 0)} + A q Tr 0)} = A q , (D.9) 

where we used the normalization condition Tr {g(t, 0)} = 1 and that the first two 
contributions on the right are null. Average values for two magnon operators are 

(c qa c q Jt) = Tr {c qo c q() 10,0)} 

= Tr { ( tqa < + « qa ^! qa + A qa ) (. qb < + Kqb # qb + A q6 ) ftt, 0)} = 

= 6 qa 6 q6 Tr {« |(t, 0)} + Vqi ,Tr {<^l q6 0)} + iqa A qb Tr {<*f qa 0)} + 

+ Vqb Tr {€ qa < 0)} + Vqb Tr {^! qa ^l q6 10, 0)} + K qa A qb Tr jVl qa |0, 0)} + 

+A qa/qb Tr 0)} + A qa K q6 Tr {<# qb £(t, 0)} + A qa A qt Tr {%{t, 0)} = 

= 6 qaKq6 Tr{<Xl q6 ^,0)}+K q/q6 Tr{^! qa < 0)) + A qa A qb Tr 0)} = 



+ 



I —F' ' F' if U qa,-q6 

1 — e i« e q 6 — 1 



(D.10) 



after using Eqs. f lD.8j) and (1D.9j) . For q a = = q we have that 



(Tq(t) = (CqC_ q |t) 



i q /t_ q 



7 + 



Kr, I — 



q^-q 



+ (Cq|t) (Cqlt) 



(D.ll) 



ic q Jt) = Tr {c qa c qt g(t, 0)} = 

Tr { + <^_ qo + A qa ) (. qb < + Kq Xl qb + A qb ) 0)} 

^q b Tr \^ a % b fa, 0)) + ^ qb Tr «^ 10, 0)) + tq< A qb Tr 



<e(*,o)} 



+ 



+« q / qb Tr {^_ qa < 0)} + < a « qi Tr {<^ qo # qb g(t, 0)} + < a A q6 Tr |^L qa §(t, 0)} + 

+A q / qt Tr |0,0)} + A qa K qb Tr £(*,())} + A qa A qb Tr {g(t,0)} = 

= 6 qa 6 q6 Tr {^ a < |0, 0)} + < a K q5 Tr {V-^Vl^ 10, 0)} + A qa A qb Tr {%{t, 0)} = 



"qal 



+ 



\K. 



q a I 



P' t 1 _p' c w qa,qb 

e * - 1 1 — e -«»a 



(D.12) 
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and, if q a = q 6 = q, 

{c\c n \t) = AT^t) = -jjsL + - K £, + |(a q |t)| 2 • (D.13) 
e i — 1 1 — e i 

Average values of more than two magnons operators are calculated in 
analogous form. 
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